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The dynamic response of asymmetric nuclear matter is studied by using a Time-Dependent Local 
Isospin Density (TDLIDA) approximation approach. Calculations are based on a local density 
energy functional derived by an Auxiliary Field Diffusion Monte Carlo (AFDMC) calculation of bulk 
nuclear matter. Three types of excited states emerge: collective states, a continuum of quasi-particle- 
quasi-hole excitations and unstable solutions. These states are analyzed and discussed for different 
values of the nuclear density p and isospin asymmetry ^ = (iV — Z)/A. An analytical expression 
of the compressibility as a function of p and ^ is derived which show explicitly an instability of the 
neutron matter around p ~ 0.09/m~'' when a small fraction of protons is added to the system. 
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I. INTRODUCTION 

The dynamic response function and the dynamic struc- 
ture factor of nuclear matter at /3-equilibrium are key 
ingredients to understand the mechanisms related to in- 
teraction processes occurring in the interior of neutron 
stars. As it is well known since the seminal works of 
Sawyer[l|,0|, the structure factor is related to the mean 
free path of neutrinos in the nuclear medium [3|. The 
opacity of nuclear matter to neutrinos, and its implica- 
tion in the neutrino emission in the early phases of exis- 
tence of a neutron star have been widely reviewed in the 
paper of Burrows et al. Q • 

Another interesting piece of information that can be 
extracted from the knowledge of the dynamic structure 
factor concerns the mechanical instability that can be re- 
lated to divergencies in the response function, and in a 
nonphysical negative value of the compressibility. This 
is particularly interesting at values of the density where 
homogeneous matter is supposed to give way to more ex- 
otic inhomogeneous phases, before reachinig the regime in 
which neutron rich nuclei are dominant [J, |5| . An inter- 
esting attempt to compute the structure factor of neu- 
tron matter by means of semiclassical simulations in this 
regime was recently performed by Horowitz et al. [6!, |7| . 

In general, several calculations of the low momen- 
tum structure factor of neutron matter have been per- 
formed using different methods, for instance in the frame- 
work of BHF-RPA theories, both at zero and finite 
temperature[8|, or in the Landau Theory approach start- 
ing from CBF effective interactions Q. 

In this paper we propose a different route to the com- 
putation of the dynamic structure factor in the isoscalar 
and isovector channels for an infinite nuclear matter with 
an arbitrary asymmetry. 
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The approach is based on first deriving a reason- 
able isospin-density functional from an equation of state 
(EoS) computed by means of the Auxiliary Field Diffu- 
sion Monte Carlo (AFDMC) ^la] method by Gandolfi et 
al.[ll|. In that paper the EoS is derived from a density- 
dependent interaction effectively in the including many- 
body contributions to the interaction. This Hamiltonian 
was fitted to correctly reproduce the saturation proper- 
ties of symmetric nuclear matter. However, the predic- 
tions of the EoS of pure neutron matter and matter at 
/3-equilibrium, including both electrons and muons, pro- 
duce realistic values of quantities such as the mass/radius 
relation in the same range of accuracy of i.e. the APR 
EoSdil. 

Obviously the range of densities that can be safely 
investigated relying on this EoS are limited. In the 
low density limit nuclear matter becomes unstable, and 
neutron matter ^hould be described including pair- 
ing correlations il,|. At densities of order 2po, where 
po=0.16fm^'^ is the nuclear matter saturation density, 
one can expect the onset of hypernuclear degrees of 
freedom [1J|, which completely changes the scenario. 

The isospin-density functional is then used in a gen- 
eralized linear response framework based on the Time 
Dependent Local Density approximation (TDLDA) . Fol- 
lowing the idea of the Local Spin Density Approximation 
(LSDA, see e.g.[l5[), and its extension to the study of the 
dynamics of the system[l6j, we derive the response func- 
tion in a Time Dependent Local Isospin Approximation 
(TDLIDA), and compute all the relevant physical quan- 
tities such as the frequency and strength of the collective 
modes and the compressibility for various values of the 
density and of the asymmetry at T=0. 

The paper is organized as follows. In Section II we will 
introduce in detail the derivation of the energy density 
functional; in Section II the extension of the TDLDA to 
the case of excitation in the isovector channel is discussed; 
Section IV contains the description of the resuts, and 
Section IV is devoted to conclusions. 
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by exactly treating the kinetic energy functional Tq. The 
detailed derivation can be found in Ref.ilSj]. The result- 
ing set of coupled self-consistent equations for neutron 
and protons wave functions, assuming ?i = c = 1, are 
given by: 



TABLE I: Coefficients of the parameterization of the 
interaction-correlation energy of symmetric nuclear matter 
(g=0) and pure neutron matter (g=l) derived by fitting the 
equation of state of Ref. [lj| . 



II. DERIVATION OF THE ENERGY DENSITY 
FUNCTIONAL 



2m 



^(r) -I- w{r)rir 



^[(r)=e,,.<^[(r), (4) 



where i stands for the set of quantunr numbers, excluding 
isospin, that characterize the single-particle wave func- 
tions. The nuclear neutron and proton densities are given 
by 



By following the Kohn-Sham methodjlTj. we introduce 
the local isospin density approximation (LIDA) for nu- 
clear matter defining the energy functional as: 



E{p,O^To{p,0+ / ev{p,0 Pdr 



(1) 



where Tg {p, ^) is the kinetic energy of the non inter- 
acting system with density p(r) — /0„(r) -I- Pp(r), and 
isospin polarization ^ = pi/p with pi{r) = Pni^) — Pp(j), 
and where p„ and pp are the neutron and proton den- 
sities, respectively. Note that Tq is the kinetic energy 
of the non-interacting system, therefore missing the con- 
tribtions coming from dynamical quantum correlations. 
However, it has a comparable magnitude and is consis- 
tently treated in this method. Equation ^ defines the 
interaction-correlation energy per particle ey(p, £,) of the 
asymmetric nuclear matter with density p and isospin 
polarization ^. This quantity can be extracted by any 
independent calculation of the total energy per particle 
of asymmetric nuclear matter after subtraction of the free 
kinetic energy contributions at each value of the p and ^. 
In this work we will use for ey (p, ^) the following 
parametrization based on the Auxiliary Field Diffusion 
Monte Carlo (AFDMC) calculations of Ref. ^ 



ey (P, = eo (p) + C^ [ei (p) - eo (p)] 



(2) 



where 



(3) 
are the interaction correlation energies of symmetric nu- 
clear (^ = 0, g = 0), and neutron matter (C = l,? = I), 
respectively. We have assumed the saturation density 
value po = 0.16fm~'^. The values of the parameters in 
Eq. ([3]) we have extracted from the fit arc given in Table 
I. 

The parametrization ([2]) reproduces very well the 
AFDMC calculations in a wide range of density p (from 
po/2 to 3po) and polarizations ^. 

By minimizing the energy functional ^ with the con- 
straint that the number of neutrons and protons remains 
constant, we can derive a set of self-consistent equations 



Pr=Y. 



^I(r)|^ 



with rjr = 1(— 1) if T = n{p) and n,p stands for neu- 
trons and protons, respectively. The effective potentials 
in Eq |4] are derived from the interaction-correlation en- 
ergy functional deriving with respect to the density and 
isospin polarization: 



viv) 



dpev[p{Y),£] 
dp{r) 



w{v) 



d^v[p{r).i)] 

ae(r) 



(5) 



As it is well known, the Kohn-Sham method sketched 
here gives an exact solution of the variational principle 
which minimizes the energy functional ([T]). The theory is 
still approximate in the sense that the exact interaction- 
energy functional is unknown, and one needs to rely 
on some expression derived by independent calculations. 
Concerning our choice for ey(p, Oi ^^ want to point out 
once more that, as discussed in Refllll the functional ([1]) 
gives an equation of state of asymmetric nuclear matter 
providing realistic predictions for neutron stars proper- 
ties when the /^-equilibrium condition is imposed. We can 
therefore be confident that both the density and isospin 
polarization dependence are sufficiently accurate to yeld 
reasonable values also for the first and second derivatives 
of ev with respect to p and ^. These are the main ingredi- 
ents of the Kohn-Sham method and of its time-dependent 
version described in the next section. 

Finally, we notice that when applied to infinite nuclear 
matter, the static equations ^ are satisfied by plane- 
wave solutions for ^pJiv), since in this case all the densi- 
ties are constant as well as the density dependent poten- 
tials ([5]). As a consequence, the Kohn-Sham equilibrium 
density and energy per particle are obviously the same of 
the starting interaction used to derive the energy func- 
tional. Hence, the static Kohn-Sham equations do not 
give any information in infinite systems. However, as de- 
scribed in the next section, the time- dependent version 
of them give new and useful solutions for the nuclear re- 
sponse even starting from the same simple ingredients. 



III. TDLIDA RESPONSE OF INFINITE 
ASYMMETRIC NUCLEAR MATTER 

In this section we calculate the TDLIDA density- 
density (x'*(9,w)) and isovector-density/isovector- 
density {x'"{q,i^)) response functions of a 3-dimensional 
spin-unpolarized uniform gas of N neutrons and Z 
protons {N -\- Z = A), with isospin polarization 

£. = Pi/P- 

We start writing the time-dependent Kohn-Sham (KS) 
equations in an external, time-dependent field along the 
r-direction: 



([9]) have solutions describing density fluctuations in the 
density operator F = J2k=i e*'*'"'', and isovector-density 
operator F^ — X]fc=i e*'*""''?/^, given by 

SF{d,io) = {i;{t)\Fm))~{0\F\0} = 
dre'^''[p{r,t) - p]^ = Ve'^'{6pn + Spp) , (10) 



and, 



SFr{d,Uj) = {ij{t)\Frm)) ~ {0\Fr\Q) = 



y^ ;^'=(g«(q-rfc-'^t) _)_ g-i(q-rfc-wt)-j 



(6) 



fc=i 



with A^ = A for density excitations and A^ — Xrjr for 
isovector-density excitations. The KS equations read: 



y'dre^'i-[pi(r,i)-pi]o 



Ve'^\5pr^-5pp) , 



(11) 



+w[p,, (r, t) , pp (r, i)]7?, + A, [e^(q--t) (7) ^-{q, u;) 



where V is the volume of the gas and O — J2k=i ^ ^^'^'' ■ 
The density-density response is given byjl8| 



A 



x"(g,c.)+xP(g,c.), (12) 



In the uniform 3D nucleon gas, the nuclear density 
oscillations induced by the external field are given by: 



Pn (r, t) = pn + Spn (r, t) , 

Pp (r, t) = pp + Spp (r, t) , 



(8) 



where pn and pp are the neutron and proton constant 
densities of the unperturbed initial state, respectively, 
and 

Spp{r,t) = 5pp(e'(i-"-"*)+e-^('i-"-"*)) , (9) 

as follows from translational invariance. The quantities 
Spn and Spp are constants to be determined. Eqs. (O and 



and the isovector-density-isovector-density response is: 

;^^(,,^) = ]^felM^^«(,,^)_^P(,,^). (13) 

In order to determine 6pn and 6pp, we then insert 
p„(r,i),pp(r, i) of Eqs. ([8l)-(l9]) into ^ and linearize the 
equations. This procedure determines the self-consistent 
KS mean-field potential entering Eq. ([7]) to be: 



VKs[Pn{r,t),Pp{r,t)] = VKs{Pn,Pp) + 
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OVks 
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Spnir,t) 



(14) 
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5pp{Y,t) 



Evaentually, from Eqs. ([7]) and ([H]) we obtain: 



1^(^,7 (r, t) = j-^V^ + const. + [5p„K.n + 5ppK,p + A] (e^(q"--"*) + ^-K^-^-^')^ \ ^^ (r, t) , 
^^<^f (r,t) = j-^V^ + const. + [<5p„K.p + <5ppK,n ± A] (e'^'^-'-'-*) + e-'^i—"*)) | ^f (r,t) 



(15) 



where in the second equation one must keep the 
plus sign when calculating the density-density response 
and the minus sign when calculating the isovector- 
density/isovector-density response. In Eqs. (|III|) we have 



defined the mean-field potentials: 
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Eq. pil[) can be rewritten as 
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- 2^V2 + C + a; [e'(i-"--"*) + e-^^'i--"*)] I (p^' (r, t) 



a coupling constant A' given by Eq. (TTS]) . For such a sys- 
tem, the density response functions are the single-particle 
free responses Xg((7,a;), Xoil^^)- From Eqs. ([T2l) and 
([T^ and from the analogous relations for the free re- 
sponse functions 



4^nr,i) 
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- V^ + C + a; [e'(i--'^*) + e-»(q--— 0] 1 ^^p (r, t) 
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we obtain 



where C is a constant, and 



K = SpnVn^n + SppVn^p + A 



(18) 






(20) 



Equation (fT7|) coincides with that of a non-interacting The solution of these equations, finally gives the TDLIDA 
system coupled to an external time oscillating field, with response functions 
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{V - Vp,pX^o)iV - K.nXSl) - Vn,pXWp,nX^O 
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{V - Vp,pXf>){V - K,nX^) - K,pXS^p,«xg 
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(21) 



Eqs. (|2T|) allow for studying the response of partially 
isospin-polarized nuclear matter (TV ^ Z), which is the 
aim of the present work. 

Note that for fully isospin-unpolarized systems (N=Z) 
the above equations are drastically simplified. In fact, 
in this case, one has: p„ = Pp = p/2, Xo = Xo = Xo/2 



and Vn 

/ dv 
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P di 
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and one gets from Eqs 
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(22) 



and 



X"(«,^) = 



xo(g,w) 



1 



_1 Oil' 
P 9? 



Xo(g,' 



(23) 



p,«=o 



which have the same form of the RPA responses of isospin 
unpolarized systems, but not the same meaning, since 
they have as main ingredients the derivative of the self- 
consistent KS potentials, and not the Fourier transform 
of an effective nucleon-nucleon interaction as in the RPA 
theory. 



The same simplifications occurs in the case of pure 
neutron matter. In this case one finds that Xo = 0- Con- 
sequently: 



x^iq,^) = x"iq,^) 



Xo(9>w) 



1-14, 



•■ V 



(24) 



Since the time-dependent-density functional approach 
only holds in the low-g, low-o; limits |18j. in the following 
we will use as the free response functions Xo ^^^ Xo ^^~ 
tering Eqs. ((2T|) the following simple expressions valid in 
these limits: 



Xo (q.^) 



-V iy""P 



■In 
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2(l±C)l/3 s+(l±C)l/3) 



(25) 



where i/"'^ 



mkp'P /tt'^ 



mkpil ± 0^ A , kp 



(^p)^/'^, s — u}/{qvf). The plus sign holds for Xo ^^'^ 
the minus sign for Xo- Note that in the low-g, low-w lim- 
its, Xo 3'ii'i Xo depend on q and lo only through the com- 
bination s = u)/{qvf). Since x*'^ of Eqs. (PT|) depends on 
q and u only via the dependence of Xo E^nd Xo on these 
variables, also the interacting responses turn out to be 
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FIG. 1: Adymensional interaction parameter Gp = VpVpp as 
a function of the isospin asymmetry ^ = (A'' — Z)/A for three 
different values of the density p. 



FIG. 2: Same as in Fig. I for the interaction parameter G„ = 



functions of s only. As in Landau theory, also in TDL- 
IDA the nuclear responses are functions of s = uj/{qvf) 
and not of q and to separately. 



IV. RESULTS 



A. Mean-field potential 



From equations (j2ip . by taking the immaginary part 
of X, it is possible to calculate the excitation strength 
S''''"{q,uj) = —{l/7r)Im x^'^ and the moments ml.'^ of the 
density (s) and isovector density (v) excitation operators 
F''", with F" = X^fe^i e'''''"" and F" = X^Li 



gjqrfc fe 



The moments are given by: 



dujuj''S'''''{q,Lo) 



E' 



<I(0|F^'1«)P.(26) 



By setting the mean- field potentials Ki„, Vpp, Vnp equal 
to zero it is possible to write: 



Sfreeil'^) 



VrnkF s 



Q(i-(TW7^)+Q(i-(T3|f7.) 



(27) 



which gives the one-particle-one-hole excitation strength, 
and, by integration, the Fermi gas energy moments of the 
asymmetric non interacting nuclear matter. 

In Figs. 1-3 and Tables II- IV the values of the adimen- 
sional interaction parameters G„ = VnVnm Gp = VpVpp 
and Gnp = ^JVnVpVnp are reported as a function of the 
isovector polarization f — [N — Z)/A for three values of 
the density p. These quantities show a strong ^ and p 
dependence, giving rise to excitations strengths and mo- 
ments which are quite different from the non interacting 
cases. 
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FIG. 3: Same as in Fig. 1 for the interaction parameter Gnp 

^JUnl^p Vnp . 



B. Energy weighted sum rules 

We first consider the moments ?n_i,r?T,i, ?7i3 of the 
density (s) and isovector density (u) operators, F" and 
F'" . These quantities still have an analytical expressions, 
given by: 



, „ _ V Vn{l + Gp) + Vp{l + Gn) T 1^/VpU;,G, 
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(30) 



where the minus sign appHes to density excitations and 
the plus sign to the isovector-density ones. 

The m^i moment is known as the hydrodynamic sum 
rule because in the N = Z case it is directly related 
to ordinary sound velocity. The mi moment is the /- 
sum rule, and the 7713 moment is related to the elastic 
properties of the Fermi systems (for a wide illustration 
of the properties of these sum rules see Ref JS). In the 
case of pure neutron matter, the m_i moment is related 
to the neutron matter compressibilty if" by the relation: 



^0 



1 



1 + G„' 



(31) 



where mP_i = Vvnil and K^ — 9n^m/k'^% respectively 
are the Fermi gas static polarizability and compressibil- 
ity. Similar expressions also hold in the symmetric nu- 
clear matter (N = Z), where !/„ = I'p, and Gp — Gn — G. 
In this case the moments read 
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(32) 



(33) 



np 



for the isoscalar and isovector cases, respectively. Expres- 
sions (I31I32P are quite similar to the results of Landau 
theory for neutron and symmetric nuclear matter[l8|. 
The true novelty of Eqs. ^ and (|28|29|30p stands in 
their explicit dependence on N, Z which allows to study 
the properties of the nuclear response when a fraction of 
protons is present in the neutron matter as happens in 
the neutron star interior. 

Note that in the local isospin density approximation 
used here, the /-sum rule mi is the same for both the 
isoscalar and isovector excitation operators, whereas cal- 
culating nil directly from the original interaction used to 
deriving the local energy functional, woud have yelded 
for ml an interaction contribution. This reflects a fail- 
ure of the model in reproducing quantities in the high cj 
region. The sum rules (|28l29l30p are only valid in the 
low-q and low-w limits, where the TDLIDA is expected 
to work. They only account for one-particle-one-hole and 
collective excitations. Many-particle-many-hole excita- 
tions are important in the high energy part of the exci- 
tation spectrum [19| , and give an important contribution 
to m\. Following this argument, since the to_i moment 



is mainly determined by the low energy part of the spec- 
trum, one expects that this sum rule should be the closest 
to the the exact value that might in principle be com- 
puted by directly solving the many-body Schroedinger 
equation. 
C. Collective modes and dynamic structure factors 

We now turn to the interacting excitation strength 
S{s — uj/qvp)- The interaction produces new types of 
excitations beyond the usual one-particle-one-hole ones 
which are the only excitations of the non interacting case. 
These excitations are given by the poles of the response 
functions (pTjl . which are the solutions of the equations: 



(1 + GpnP){l + Gnfl'') - GlpQ^^QP = 



(34) 
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2(l±4)l/3 ^^^ s+(l±C)l/3) 

sign holding for il" and the minus sign for 17^. Depend- 
ing on the strength of the interaction, these new solu- 
tions are essentially of three types: a) real solutions such 
that s > (1 ± ^)^/^ (collective modes), producing a dis- 
crete peak in the dynamic form factor S{s) with no at- 
tenuation; b) Solutions with some imaginary component, 
and corresponding modes which decay by exciting single 
quasi- particle-quasi-hole pairs (Landau damping) ; c) Un- 
stable solutions which are associated to the divergence of 
the polarizability sum rules ml^^- 

In Tables II-IV, we report the values s of s at which 
the collective states occur in the isoscalar (sg ) and isovec- 
tor (si,) strengths, the percentage mf and m^" of energy- 
weighted sum rule exhausted by the collective states in 
the isoscalar and isovector channels, respectively, and the 
two mean s values 
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(35) 



The energy uJi^i is here- 



obtained by Eqs. P8I29I30|) . 
after referred to as the hydrodynamic energy in analogy 
to what happens in liquid He'^, where vi ~ ojf _i/q re- 
produces the predictions of the hydrodynamical model 
for the ordinary (first-) sound wave. Conversely the en- 
ergy Wg'J' is referred to as the elastic energy since in liquid 
He"^ vq — uj^i/q reproduces the predictions of the elastic 
model for the zero-sound wave[l8|. 

The strength of the collective state at s = s can be 
computed from the following (adymensional) expression 
for the response: 
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(36) 
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FIG. 4: Quasi-particle-quasi-hole excitation strengtiis in units 
of Nm/2kp as a function of s = u/qvp for three differ- 
ent values of the density (p = 0.30/m~^, p — 0.16/m~'^, 
p = 0.10/m~"^) at fixed isospin asymmetry ^ = 0. The posi- 
tions of the collective states, when present, are indicated by 
an arrow. The dashed and full lines stand for the isoscalar 
and isovector excitation strengths, respectively. The fraction 
of energy-weighted sum rule exhausted by the two strengths 
is also reported near the curves. The remaining fraction of 
energy-weighted sum rule is exhausted by the collective states. 



FIG. 5: Same as in Fig. 4 at isospin asymmetry ^ = 0.3. The 
collective state, indicated with an arrow, have the same en- 
ergy both in the isoscalar and isovector channels. The dashed 
and full lines stand for the isoscalar and isovector excitation 
strengths, respectively. The fraction of energy-weighted sum 
rule exhausted by the two strengths is also reported near the 
curves. The arrow is dashed or full depending on whether the 
collective state exhausts more strength in the isoscalar or in 
the isovector channel. 



derived from Eqs. ((2T|) and (|25|) by expanding x*'" 
around the pole at s. Naming N{s) and D{s) the nu- 
merator and the nenominator of Eq. (|36p . respectively, 
one obtains: 



S{s) 



Nm/{2kj,) 



dPjs) 
ds 



6{s 



(37) 



The fraction of the /-sum rule mi exhausted by the col- 
lective state is easily calculated to be 



N{s) 



aD(s) 

ds ' 



(38) 



In Figs. 4-8 we then plot the quasi-particle-quasi- 
hole excitation strengths S''*'^(s — u/qvp) in units of 
Nm/{2kp) for different values of the density p and s at 
fixed values of the isospin polarization ^. This quantity 
has been obtained by numerically computing the immagi- 
nary part of expression (|36p . In the figures the fraction of 
energy weighted sum rule exhausted by the continuum of 
quasi-single particle states is explicitly indicated, and the 
position of the collective state, if present, is indicated by 
an arrow. In all the calculations presented in the follow- 
ing, we have numerically checked that the particle-hole 



^ = 0.8 



p = 0.10fm" _ 



p = 0.16 fm" 






p = 0.30 fm " 



s = m/qv 
FIG. 6: Same as in Fig. 5 at isospin asymmetry ^ — 0.! 



and collective contributions to the strength completly ex- 
haust the mi sum rule. Let us start the discussion of the 
two extreme cases corresponding to ^ = (symmetric 
nuclear matter N = Z) and ^ = 1 (pure neutron matter 
Z = 0). In the first case {N = Z), isoscalar and isovec- 
tor modes are decoupled in the sense that the isoscalar 
and isovector density operators F* and F" give rise to 



c 


Ss 


% ml 


Sv 


%ml 


<-i 


Sl.-l 


^3,1 


^h 


G, 


G„ 


Gnp 


0.0 


- 


- 


1.1036 


81.4 


0.115 


0.915 


0.529 


1.051 


0.2757 


0.2757 


-1.2363 


0.3 


1.1154 


7.7 


1.1154 


58.2 


0.142 


0.925 


0.472 


1.088 


0.3991 


0.1422 


-1.2031 


0.5 


1.1466 


5.5 


1.1466 


15.1 


0.183 


0.890 


0.472 


1.131 


0.4611 


0.0561 


-1.1409 


0.8 


- 


- 


- 


- 


0.269 


0.891 


0.564 


1.105 


0.4690 


-0.0599 


-0.9578 


0.9 


- 


- 


- 


- 


0.313 


0.915 


0.632 


1.129 


0.4093 


-0.0937 


-0.8408 


1.0 


- 


- 


- 


- 


0.681 


0.681 


0.942 


0.942 


0.0000 


-0.1247 


0.0000 



TABLE II: Numerical results at density p — O.lOfm^"^ and various isospin polarizations for: a) the energy s of the collective 
mode in the isoscalar (s) or isovector (v) channels, and corresponding fraction of the /-sum rule exhausted by such modes; 
b) the centroids of the collective peaks excitations predicted by the energy- weighted sum rules (see Eq. (34)); c) Interaction 
parameters giving the mean-field potential. 
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Ss 


%mf 


Sv 


%m1 


sf.-i 


^i.-i 


4,1 


4,1 


G, 


G„ 


Gnp 


0.0 


- 


- 


1.1008 


80.9 


0.609 


0.911 


0.799 


1.047 


0.8012 


0.8012 


-0.6877 


0.3 


1.1316 


17.8 


1.1316 


60.7 


0.614 


0.909 


0.793 


1.076 


0.8759 


0.7077 


-0.6610 


0.5 


1.1685 


26.8 


1.1685 


46.9 


0.623 


0.904 


0.820 


1.102 


0.8953 


0.6464 


-0.6125 


0.8 


1.2281 


29.6 


1.2281 


33.6 


0.650 


0.901 


0.916 


1.150 


0.7997 


0.5669 


-0.4826 


0.9 


1.2481 


29.1 


1.2481 


30.7 


0.667 


0.908 


0.967 


1.167 


0.6751 


0.5455 


-0.4106 


1.0 


1.2680 - 


28.5 


1.2680 


28.5 


0.899 


0.899 


1.110 


1.110 


0.0000 


0.5271 


0.0000 



TABLE III: Same as table II for density p = 0.16fm" 




s - CO/qv 
FIG. 7: Same as in Fig. 5 at isospin asymmetry ^ — 0.9. 



distinct isoscalar and isovector strengths S"'(s — uj/qvp) 
and S^{s — uj/qvp) ■ This is different from what happens 
when N y^ Z where F* and F" can indifferently excite 
modes both in the isoscalar and in the isovector channels, 
with strengths and weights depending on the values of ^ 
and p. Obviously, for neutron matter (Z = 0) it makes 
no sense to speak of isoscalar and isovector modes since 
in this case there is only one type of excitation. 

For symmetric nuclear matter (^ = 0) at density 
p = O.SOfm^'^, in the isoscalar strength S'^{s) a strong 
collective state is present at Sg = 1.1188, exhausting 
about 83.6% of the energy weighted sum rule mi. A 




p = 0.10fm" 

^=1.0 



p = 0.16 fm" 



p = 0.30 fm 



■- m/qVp 



FIG. 8: Same as in Fig. 5 at isospin asymmetry ^ = 1. 
In this case, pure neutron matter, there is only one type of 
excitation and the distinction between isoscalar and isovector 
channel does not apply. 



continuum of single particle type excitations is instead 
predicted at lower s values. A similar situation occurs 
for S'"{s). Here the collective state is at at Sy = 1.0958, 
and exhausts about 80.1% of the energy weighted sum 
rule TOi. The collective states in both the cases occurr at 
energies which are closer to the mean excitation energies 



^3,1 



than to the 



ones, showing that the collective 



modes are of mainly of elastic type. This situation is 
analogous to what happens, for instance, in liquid He'^. 
By decreasing the density, we still predict the existence 
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Ss 


%ml 


Sv 


%m1 


Sl.-l 


Sl.-l 


^3,1 


^3.1 


Gp 


G„ 


Gnp 


0.0 


1.1188 


83.6 


1.0958 


80.1 


0.937 


0.903 


1.070 


1.040 


1.5414 


1.5414 


0.0935 


0.3 


1.2072 


58.2 


1.2072 


48.0 


0.943 


0.902 


1.099 


1.059 


1.5260 


1.5267 


0.1091 


0.5 


1.2643 


64.3 


1.2643 


57.8 


0.953 


0.900 


1.148 


1.094 


1.4732 


1.5146 


0.1351 


0.8 


1.3420 


74.6 


1.3420 


71.7 


0.986 


0.902 


1.263 


1.188 


1.2229 


1.5064 


0.1829 


0.9 


1.3666 


78.0 


1.3666 


76.5 


1.011 


0.909 


1.310 


1.232 


1.0106 


1.5084 


0.1893 


1.0 


1.3907 


81.4 


1.3907 


81.4 


1.153 


1.153 


1.324 


1.324 


0.0000 


1.5134 


0.0000 



TABLE IV: Same as table II for density p = 0.30fm" 



of a strong collective state in the isovector channel, ex- 
hausting about 80% of the energy weighted sum rule mi. 
On the other hand, the collective mode disppears in the 
isoscalar channel, due to Landau damping. 

The situation is similar, but not exactly the same, in 
the case of pure neutron matter (^ = 1). At high den- 
sity (p — 0.3 fm~'^) the dynamic form factor shows a 
strong collective state of elastic nature at s = 1.3907, ex- 
hausting about 81.4% of the energy weighted sum rule, 
together with a continuum of single particle excitations 
at lower s values. Decreasing the density to p = 0.16 
fm~'^, the energy of the collective state decreases, like- 
wise the percentage of rrii sum rule exhausted by this 
state. Eventually, at some density below p — 0.16 fm^'^, 
the collective completely decays in one particle-one hole 
excitations. 

When N ^ Z, isoscalar and isovector probes excite an 
admixture of isoscalar and isovector modes. In all cases 
we predict a single collective state, when present, having 
different weights in mi for the isovector and isoscalar 
channels. 

The typical situation expected in the neutron star core 
is that of a small quantity of protons present in the neu- 
tron matter, corresponding to ^ close to 1. Rather than 
considering the proton fraction predicted for each value 
of the density by the AFDMC equation of state, we have 
considered the cases £, — 0.8 and £, = 0.9 for all den- 
sities, in order to have an idea of the evolution of the 
collective modes with this parameter. For such values 
of the asymmetry, a clean collective state in both chan- 
nels is present only at high density, where it exhausts 
about 70-80% of the rrii sum rule. Already at saturation 
density p — 0.16fm~'^, this state has almost completly 
decayed into particle-hole excitations, since it exhausts 
only about the 30% of the rrii sum rule in both channels. 
At p = 0.10/m^'^ it is completly damped. This sug- 
gests that collective modes in the isoscalar or isovectors 
channels can be present going inward from the interface 
between the inner crust and the outer core. 

Another interesting analysis can be made by looking at 
the compressibility (in unit of the free one) as a function 
of the density for different values of the proton fraction. 
This quantity is given by m'Li/m^_i and m'"_i/'mP_i, and 
is plotted in Figs. 9. It can be noticed that a strong 
divergence for a density near to p = 0.085fm~^ appears at 
^ = 0.9, and that it moves to a larger density increasing 



the proton fraction. We interpret this divergence as a sort 
of mechanical instability of matter towards the formation 
of an inhomogeneous phase, as expected in the inner crust 
of the neutron star. 



It is also interesting to look at the regime that would 
correspond to asymmetries typical of very neutron rich 
nuclei (^ = 0.3). At (unphysical) high density, a col- 
lective state exhausting about the 50-60% of toi at 
s = 1.2072 is present both in S''{s) and S"'(s). At the 
same time the quasi-particle-quasi-hole strength is prac- 
tically concentrated in one state at a value of s sligthly 
smaller than 1 and exhausting about the 50-40% of mi . 
By lowering the density, we observe that at p = 0.16fm~'^ 
the collective state appears at a lower energy, and re- 
mains collective only in the isovector channel, where it 
still exhausta more than 60% of theenergy weighted sum 
rule. In the isoscalar channel most of the strength is 
taken by the one-particle-one excitations, though con- 
centrated in a narrow region of s, and the "collective" 
solution exhausts only about 17.8% of the mi sum rule. 
At p = O.lOfm"'^ the collective mode can be excited only 
by an isovector probe, and in the isoscalar channel it is 
practically completly damped. The quasi-particle-quasi- 
hole strength is distributed in a wide range of s in the 
isoscalar channel, whereas in the isovector one continues 
to see a substantial concentration at s ~ 1. 



This situation might be interpreted in analogy to what 
happens in nuclei with a large excess of neutrons, where 
a low energy peak beyond the usual giant resonance is 
observed in the isovector channel. 



Finally, in Figs. 10 we plot the compressibility ratio for 
three different values of the density. At ^ = and ^ = 1 , 
m''_i/'mP_i gives the compressibility of symmetric nuclear 
matter and pure neutron matter, respectively. These val- 
ues are not far from the results computed by other au- 
thors starting from a microscopic AV8' Hamiltonianj2(| . 
From the figures one also sees that differently from what 
happens at low density, at normal and high densities, 
ra'^__i/mP_i is practically independent on ^. On the con- 
trary, m'"_i/'wP_i results to be always an increasing func- 
tion of the isospin asymmetry. 
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FIG. 9: Isoscalar and isovector compressibility ratios 
mLi/m^Li and m'Li/m^i as a function of the density p for 
values of ^ close to 1. 
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FIG. 10: Isoscalar and isovector compressibility ratios 
Tnli/jn'li and m!Li/m'li as a function of the isospin asym- 
metry 5 for three different values of the density. 



V. CONCLUSIONS 

We have studied the dynamic form factor of asym- 
metric nuclear matter by using a time-dependent local 
isospin density approximation approach based on a local 
density energy functional derived by an Auxiliary Field 
Diffusion Monte Carlo calculation. The more relevant 
results we have found are the foUovifing: i) The presence 
of a strong collective state at high density at an energy 
which increases with the value of the nuclear asymmetry 
^ — {N — Z)/A; decreasing the density down to sat- 
uration {p = 0.16 fm~^), this state tends to decay in 
quasi-particle-quasi hole states. The decay is faster in the 
isoscalar channel, where it remains very collective only at 
values of ^ close to 0. Further decreasing the density, the 
collective states survives only in the isovector channel at 
small values of ^, which are typical of the neutron star 
interior, ii) When a small fraction of protons is added 
to the neutron matter and ^ is equal to 0.8, 0.9, at val- 
ues of the density near or slightly smaller of 0.09 fm~^ 
the system becomes unstable. This unstability is seen 
in the isoscalar and isovector compressibilities, which at 
such densities diverge, iii) When ^ is around 0.3 (small 
asymmetry), in the isovector channel at all the densities 
two states cohexist, one of collective and the other of 
quasi-particle-quasi-hole nature, practically sharing the 
fraction of exhausted energy-weighted sum rule. This is 
similar to what is observed in the photodisintegration of 
large neutron excess nuclei. 
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